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We present initial results of the use of Force Gradient integrators for lattice field theories. These 
promise to give significant performance improvements, especially for light fermions and large 
lattices. Our results show that this is indeed the case, indicating a speed-up of more than a factor 
of two, which is expected to increase as the integration step size becomes smaller for larger lattices 
and smaller fermion masses. 
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1. Introduction 

Our goal is to construct more accurate molecular dynamics integrators for use in Hybrid Monte 
Carlo (HMC) computations for lattice quantum field theories. In particular we shall consider sym- 
metric symplectic integrators whose errors are of higher order in the integration step size than those 
of the leapfrog (also known as the Stormer or Verlet) method. These integrators improve the scaling 
behaviour SH = k ■ 8x" from n = 2 to n = 4, which reduces the cost for large enough volume and 
small enough fermion masses for which 8t — > 0. Previous methods for constructing higher-order 
integrators ^ |2|] were thwarted by large values for the coefficient k; the new idea considered here 
is to compute second derivatives "analytically" rather than "numerically" [|], Q]. 

2. Symplectic Integrators 

As described in [gg^Jl we may define a Hamiltonian system for a gauge field by introducing 
the symplectic fundamental 2-form co = —d(p'di) with 0, being the frame of left-invariant Maurer- 
Cartan forms; this ensures that the Hamiltonian dynamics is gauge invariant. For every 0-form F 
on phase space this defines a Hamiltonian vector field F satisfying dF = ipCO, and the Hamiltonian 
evolution for the system corresponds to an integral curve of the Hamiltonian vector field H for the 
Hamiltonian function H. We can find a closed-form integral curve of F, that is evaluate e F T explic- 
itly, if F depends only on the "positions" (fields) q or momenta p. This is particularly useful when 
the Hamiltonian is of the form H(q,p) = S(q) + T(p) as then we can integrate the Hamiltonian 
vector fields 5 and f exactly. 

3. Shadow Hamiltonians and Force Gradient Integrators 

We recall the Baker-Campbell-Hausdorff (BCH) formula, which states that if A and B belong 
to any (in general non-commutative) associative algebra then e A e B = e A+B+s where 8 is in the 
free Lie algebra generated by A and B. Furthermore it gives an explicit expansion for 8, and 
correspondingly for the symmetric product 



Using the Jacobi identity one may show that the commutator of two Hamiltonian vector fields 
is itself a Hamiltonian vector field, [S, f] = {S, T}, where {S, T} = —co(S, f) is the Poisson bracket 
of the two 0-forms S and T. We therefore find that any integrator constructed from a sequence of 
symplectic steps exactly conserves a shadow Hamiltonian H obtained from the BCH formula by 
replacing commutators with Poisson brackets. 

As a very simple example consider the PQPQP integrator 




(\A,\A,B]]+2[B,[A,B}]) + 




whose shadow Hamiltonian is 



H = H + 



6a 2 - 6a + 1 



{Ms,r}} + 



1 -6a 



{T,{S,T}} )8t 2 + 0(8t 4 ). 
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Shadow Hamiltonian 
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T + S - 6 4 ({S, {S, T}} + 2{T, {S, T}})+ @{8x A ) 


QPQ 




T + S+ s 4 (2{5, {S, T}} + {T, {S, T}})+ ^(St 4 ) 


PQPQP 

a=l 

§§§ 


e \bxS e \8xf 


T + S+^{S,{S,T}} + 0(5t 4 ) 


PQPQP 

§§§ 


e ^l8xS e \Sxt 
x e^ 
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Alt 



x e 



r+s 

/ -(40^4+40^2+48) {5, {S, {S, {S, T}}}} \ 
+ (180^4+240^2+312) {{S, T},{S, {S, T}}} 
+ (60^4+80^2+104) {{S, T}, {T, {S, T}}} 

+ (-20^4+8) {T, {S, {S, {S, T}}}} 

+ (20^2+32) {T,{T, {S, {S, T}}}} 

\ +(5^2+8) {T,{T,{T,{S,T}}}} > 
+^(5t 6 ) 
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Force 
Gradient 
PQPQP 



48<StS-gT 3 {S+S/T}} 

x e 72 



T + S 

( 41 {S, {S, {S, {S, T}}}} \ 
+36 {{S,T},{S,{S,T}}} 
+12{{S,T},{T,{S,T}}} 
+84{T,{S,{S,{S,T}}}} 
+ 126 {T,{T,{S,{S,T}}}} 

\ +54 {T,{T,{T,{S,T}}}} ) 



+ 



Table 1: A selection of integrators with their exactly conserved shadow Hamiltonians. 



As the Poisson bracket {S, {S, T}} does not depend on momentum we can integrate the Hamil- 
tonian vector field {5,{5,r}} exactly, and this "second derivative" corresponds to the Force Gra- 
dient just as the Hamiltonian vector field S corresponds to the "force". The explicit form of the 
shadow Hamiltonian for a variety of integrators is shown in Table [j], the simplest Force Gradient 
integrator is given in the last entry. 



4. Computing the Force Gradient 

As shown in our previous proceedings ||] Poisson brackets for gauge theories may be writ- 
ten in terms of momenta pj and linear differential operators e, that provide gauge-covariant gener- 
alizations of the vector fields d/dq l ; in particular we have {S, {S, T}} = e'(S)ei(S). The "equations 
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log 10 6t 

Figure 1: Change 8H in the Hamiltonian over an entire trajectory as a function of the integration step size 
8t. The dashed line indicates the value of 8H which corresponds to a 75% HMC acceptance rate. 

of motion" for the {S, {S, T}} vector field are P = {S, {S,T}}P = —G where P = piT', and we may 
use the relation e,-(t/) = —Till, where 7} are the adjoint generators of the gauge group, to show that 
the Force Gradient vector field is the "second derivative" G = e J (S)ejei(S)T' . 

If we consider a generic pseudofermion action S = <p^^~ l (U)(j> where <j> is a pseudofermion 
field, U the gauge field, and any hermitian fermion kernel, then 

ei (S) = -tf Jt~ X ei{J()J(-^ = -X*ei{JZ)X 

where X = The Force Gradient can be computed by applying the linear differential operator 

& = e-*(S)ej to the above equation; by the Leibnitz rule 

-G=-&[ ei {S)) = ^{X^)e i {^)X+X^{e i {^())X + X^e i {^)^{X)\ 

defining Y = &(X) = -Jt- x &{Jt)X, we obtain 

—G = Y^ ei {J()X +X t ,^(e i (^))X + X Jf e i (^)Y. 

Note that the cost of computing the Force Gradient in an HMC integrator is the inversion required 
to compute Y in addition to the usual inversion needed to compute X. 

5. Results 

We have implemented this PQPQP Force Gradient integrator for lattice QCD with dynamical 
Wilson fermions, and we present our initial results for a 16 4 lattice at j8 =5.5 and K = 0.1575, 
which corresponds to a pion mass of m K = 665 MeV. We used the usual even-odd preconditioning, 
which does not introduce any significant complications into our formalism. 
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Figure 2: Change 8H in the Hamiltonian over an entire trajectory as a function of an estimate of the com- 
putational cost. The dashed line again indicates the value of 8H corresponding to a 75% HMC acceptance 
rate. 

In Figure [l] we show the change 8H in the Hamiltonian over an entire trajectory as a function 
of the integration step size 5t on a log-log plot. To guide the eye we have also drawn a line to 
indicate the value of 8H which corresponds to a 75% HMC acceptance rate. The slopes of the 
lines correspond to the expected order of the integrators, and the Force Gradient integrator is more 
than order of magnitude more accurate than the Campostrini integrator at any step size, indicating 
that the coefficient k discussed in the introduction is indeed much smaller. 

In Figure ^ we replot the same data as a function of an estimate of the cost, namely the number 
of CG solutions divided by the step size. The Omelyan integrator requires two inversions of the 
Wilson-Dirac operator per step, one 1 for each 5 (P) integration step, whereas the Campostrini and 
Force Gradient integrators require three inversions (one for each S and {S, {S, T}}). From the 
intercepts with the dashed line (75% acceptance) we find that the Force Gradient integrator is a 
factor of 2.6 cheaper than the Omelyan integrator, which was up to now the preferred choice of 
integrator. 

6. Conclusions 

Our Force Gradient integrator is cheaper by more than a factor of two even for small lattices 
with fairly heavy quarks, and the benefit increases as the integration step size becomes smaller. 
We expect that the integrators will be improved by tuning using measured average values of Pois- 
son brackets, as described in We also note that our formalism for computing Force Gradient 
integrators is compatible with all common actions, smearing, and so forth. 

1 As initial and final "half steps" can be combined we count each as half an inversion. 
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